Ah initio electronic structure calculation of hoUandite vanadate K2V8O16 
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An ab initio electronic structure calculation based on the generalized gradient approximation in 
the density functional theory is carried out to study the basic electronic states of hoUandite vanadate 
K2V8O16. We find that the states near the Fermi energy consist predominantly of the three t2g- 
orbital components and the hybridization with oxygen 2p orbitals is small. The dy^ and d^^ orbitals 
are exactly degenerate and are lifted from the dxy orbital. The calculated band dispersion and Fermi 
surface indicate that the system is not purely one-dimensional but the coupling between the VO 
double chains is important. Comparison with available experimental data suggests the importance 
of electron correlations in this system. 

PACS numbers: 71.30.+h, 71.20.Be, 71.28.+d 



I. INTRODUCTION 

The metal-insulator phase transition associated with 
charge and orbital ordering has been one of the central 
issues in physics of strongly correlated electron systems. 
Recently, Isobe et al} reported that the metal-insulator 
transition occurs in hoUandite vanadate K2V80ig at 
^160 K, which is accompanied by a rapid reduction of the 
magnetic susceptibility. A characteristic superlattice of 
\/2a X \/2a x 2c is observed below the transition tem- 
perature. A possible charge-ordering phase transition 
accompanied by the spin-singlet formation has thereby 
been suggestedi^ Quite recently, the phase diagram under 
high pressures has also been obtainedj^- where a variety 
of charge-ordered phases have been suggested to appear. 
The metal-insulator transition has also been observed in 
Bi.VgOie^^ and KaCrgOie-^ 

The crystal structure of K2 VsOig belongs to a group of 
hoUandite-type phases and has a VgOie framework com- 
posed of double strings of edge-shared VOe octahedra. 
The system may be regarded as a one-dimensional ver- 
sion of LiV02 known as a possible orbital-ordering sys- 
tem of the t2g orbitals on the triangular lattice of S* = 1 
spins.*' The present system K2V8O16 has the average va- 
lence of V^'^^+ and thus is in the mixed valent state of 
V3+. = idP : = 1 : 3. Thus, the central issue 
in the present system is the mechanism of the metal- 
insulator transition concerning how the highly frustrated 
spin, charge, and orbital degrees of freedom at high tem- 
peratures are relaxed by lowering temperatures and what 
type of order is realized in the ground state. 

In our previous paper^i we studied this material in 
the strong correlation limit; we set up a high-energy 
model Hamiltonian and applied a strong-coupling per- 
turbation theory to obtain the low-energy effective spin- 
orbit Hamiltonian. The obtained effective Hamiltonian 
was analyzed numerically and the possible orbital and 
spin structure of the ground state of the system was 
proposed J 

However, so far not much is known even for the basic 
electronic states of the system, such as the band disper- 



sion, Fermi surface, density of states, etc., both experi- 
mentally and theoretically. The purpose of the present 
paper is therefore to present an ah initio electronic struc- 
ture calculation of this system to clarify its basic elec- 
tronic structure. Here, we use the generalized gradient 
approximation (GGA) for the electron correlations in the 
density functional theory. The study is thus in the weak 
correlation limit but we hope that one can learn much of 
the basic electronic states of this material. 

In this paper, we will show the following: The elec- 
tronic states near the Fermi level are dominated by the 
three i2g-orbital components, where the hybridization 
with the oxygen 2p orbitals is small. The dyz and d^x 
orbitals are exactly degenerate and are lifted from the 
dxy orbital (hereafter we use the notation of the axes x, 
y, and z defined in Ref.,^ also see Fig. 5 below). The cal- 
culated band structure and Fermi surface indicate that 
the system is not purely one-dimensional but the cou- 
pling between the VO double chains is important. We 
find no nesting features in the Fermi surfaces that con- 
tribute to the instability of the doubling of the unit cell 
along the c axis, suggesting the observed lattice instabil- 
ity to be the strong-coupling origin. We will compare our 
calculated results with experiment such as photoemission 
spectroscopy^ and nuclear magnetic resonance (NMR);^ 
from which we will discuss the implications on the elec- 
tronic states and peculiarity of this material. 

We hope that the present study, if combined with the 
strong-coupling approach, will give a reliable theory to 
explain the mechanism of the electronic phase transition 
of this intriguing material. 



II. METHOD OF CALCULATION 

We here employ the computer code WIEN2kiS based 
on the full-potential linearized augmented-plane-wave 
(FLAPW) method. We have tested both the local density 
approximation (LDA) and generalized gradient approxi- 
mation (GGA) for the exchange-correlation potentia l^ 
in the density functional theory. In this paper, we will 
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however present only the results of GGA because no sig- 
nificant differences are found between the results of LDA 
and GGA. The spin-orbit interaction is not taken into 
account. No spin-polarization is assumed. In the self- 
consistent calculations, we use 4,335 k points in the ir- 
reducible part of the Brillouin zone (see Fig. 1) with an 
anisotropic sampling in order to achieve better conver- 
gence. We use the plane- wave cutoff of Xmax = 4.24 
Bohr-i. 

We assume the experimental crystal structure of 
K2V8O16 observed at room temperature (above the 
metal-insulator transition) with the lattice constants of 
a = 9.963 and c = 2.916 AM The Bravais lattice is body- 
centered tetragonal (see Fig. 1) and the primitive unit cell 
contains four V ions, one K ion, and eight O ions. We 
use the code XCrySDenii for graphical purposes. 




FIG. 1: (Color online) Schematic representation of (a) 
the unit cell of the body-centered tetragonal lattice (solid 
lines) and (b) Brillouin zone of K2V8O16. In (a), the 
primitive unit cell is also shown in the thin dotted lines. 
In (b), the symbols represent r(0,0,0), M(27r/a, 0, 0), 
X{n/a,Tv/a,0), Y'{-k / a,-K / a,-K / c), Ki(0, 0, 7r(l/c -I- c/a^)), and 
K2(27r/a, 0, 7r(l/c — c/a^)), where Ki and K2 are equivalent. 



III. RESULTS OF CALCULATION 

A. Band dispersion 

The calculated band dispersion near the Fermi energy 
is shown in Fig. 2. We find that there are 12 t2g bands 
of the V 3(i orbitals, 5 of which cross the Fermi level. We 
find the highly dispersive bands along the F-Ki, X-P, and 
M-K2 lines and weakly dispersive bands along the F-X, 
P-Ki, and K1-K2 lines, reflecting the one-dimensionality 
of the electronic state. However, the band structure is 
not simple even near the Fermi energy, suggesting that 
the contributions from the three-dimensionality are not 
negligible; in particular, the relatively large dispersion 
along the F-M line reflects the coupling between the VO 
double chains. Thus, to consider the low-energy elec- 
tronic properties of this system, the three-dimensionality 
or the coupling between the VO double chains is essen- 
tial. This result is consistent with the observed rather 



small anisotropy of the electric resistivity of the single 
crystal of this materiali^ 
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FIG. 2; Calculated band dispersion of K2V8O16 near the 
Fermi level (horizontal line). There are 12 t2g bands of the 
V 3ci orbitals, 5 of which cross the Fermi level. The labels of 
the fc-points are shown in Fig. 1. 



B. Density of states 

The calculated density of states is shown in Fig. 3, 
both in a wide energy range in (a) and near the Fermi 
energy in (b). We find that the V id band is located 
between —0.5 eV and 4.5 eV and is well separated from 
the O 2p band located between —7.5 eV and —2.5 eV. The 
hybridization between V 3d and O 2p orbitals seems to 
be rather small, justifying the use of the strong-coupling 
model with only the d orbitalsii We also find that the 
t2g band is well separated from the Cg band in the V id 
bands, corresponding to a large value of the crystal field 
splitting of lODq ~ 2 eV. Thus, the low-energy properties 
of this system is essentially governed by the three t2g 
orbitals, d^y, dyz, and dzx- 

The calculated orbital-decomposed partial densities of 
states Pa{s) {a = xy,yz,zx) are shown in Fig. 3(b). We 
first find that the total density of states at the Fermi 
level is very high, 30.4 states/f.u./eV (see Fig. 3(a)), 
which comes mainly from the three t2g orbitals; the large 
ferromagnetic spin fluctuations observed in the NMR 
experimen1>^ may be related to this high density of states 
at the Fermi level. We then find that the three t2g or- 
bitals equally contribute to the low-energy states near 
the Fermi energy. More precisely, we note that, even in 
the presence of the distortions of the VOg octahedra, the 
relation Pyz{£) — Pzx{£) strictly holds. This is due to 
the inversion symmetry along the c-axis preserved even 
in the presence of the distortionii We find that the width 
of the dxy band is slightly smaller than the width of the 
dyz and dzx bands. We also find that the number of elec- 
trons in the dxy orbital is slightly larger than that in the 
dyz and dzx orbitals; in the V atomic sphere, they are 
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0.345 and 0.246, respectively, in the energy range of the 
bands (see Sec. Ill D). 




Energy (eV) 

FIG. 3: (Color online) (a) Calculated density of states (per 
formula unit (f.u.)) of K2VgOi6 in a wide energy range, (b) 
Calculated orbital-decomposed partial density of states (per 
atom) near the Fermi level. The Fermi level is indicated by 
the vertical line. Contributions from the dyz and d^x orbitals 
(thin dotted lines) are exactly degenerate. 



We should point out that the possible nesting features 
in the two nearly-parallel Fermi surfaces seen in Fig. 4(c) 
does not lead to the 2kF instability corresponding to the 
observed^ doubling of the unit cell along the c axis ex- 
pected in the weak-coupling theory because the distance 
between the two Fermi surfaces is too short. No other 
Fermi surfaces are seen to contribute to the instability of 
the \/2a x V^o, x 2c superlattice formation. We there- 
fore consider the observed lattice instability to be the 
strong-coupling origin, such as the long-range Coulomb 
interactions leading to charge orderingsZ 




FIG. 4: (Color online) Calculated Fermi surfaces of K2V8O16. 
The 53rd to 57th bands counted from the lowest are shown 
in (a) to (e), respectively. 



C. Fermi surface 

The calculated results for the Fermi surface are shown 
in Fig. 4. There are 12 t2g bands, 5 of which cross the 
Fermi level and form the Fermi surface. If the system 
were strictly one-dimensional, a pair of the parallel Fermi 
surfaces should appear at — ifep. We actually find 
such types of the pairs of the Fermi surfaces in Figs. 4 (c) 
and (d), but the surfaces are rather distorted and more- 
over there are other types of the Fermi surfaces present. 
In particular, the Fermi surface shown in Fig. 4(e) comes 
from the dispersive band along the F-M line of the Bril- 
louin zone (see Fig. 2) and is related to the couphng 
between the VO double chains. Therefore, we may con- 
clude that the system is not purely one-dimensional, but 
the three-dimensionality or the coupling between the VO 
double chains is important. 



D. Valence electron distributions 

The spatial distribution of the valence electrons of 
K2V8O16 in the energy range between —0.5 eV and the 
Fermi level is shown in Fig. 5. The electrons mainly 
contribute to this distribution. We find that the result 
can be understood if we assume that the valence elec- 
trons are confined predominantly in the ligand field of 
each VOe octahedron; i.e., the distribution can be pro- 
jected mainly onto the d^y, dyz, and dzx orbitals. More 
precisely, however, we find that the distribution is some- 
what compressed in the z direction (toward the apical O) 
and spread in the xy plane, indicating that the number 
of dxy electrons is larger than that of dyz and dzx elec- 
trons. This is consistent with the results for the orbital- 
decomposed partial density of states (see Fig. 3), where 
the area below the Fermi level in the dxy component is 
slightly larger than that in the dyz and dzx components. 
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We also note that the electron distribution reflects the 
symmetry that the dy^ and d^x orbitals are equivalent. 
In the atomic limit, this degeneracy corresponds to the 
fact that the atomic energy levels of the dyz and dzx or- 
bitals are exactly degenerate in the ionic model and are 
located slightly lower in energy than the atomic level of 
the dxy orbitalji The orbital ordering in this system has 
thereby been suggested in the strong coupling theoryi^ 




FIG. 5: (Color online) Calculated spatial density distribu- 
tions of the valence electrons of K2VgOi6 in the energy range 
between —0.5 eV and the Fermi level, where the t2g electrons 
mainly contribute. The isovalue of 0.08 electrons/Bohr"^ is 
used. 



to higher binding energies, reflecting the opening of the 
quasiparticle gap of 230 meV<^ The reconstruction of the 
electronic states by the metal-insulator transition should 
be pursued theoretically. However, the LDA and GGA 
calculations cannot reproduce the charge-ordered states 
with the spin-singlet formation. One needs to invent the 
strongly-correlated electron models to account for the 
phase transition, of which the simplest one-dimensional 
version has been proposed in Refii. Further study is re- 
quired to take into account of the three-dimensionality of 
the system. 

It has been reported in a recent NMR experiment^ 
that there appears an exceptionally strong temperature 
dependence of the rotation of the NMR Knight-shift ten- 
sor, which suggests the strong orbital dependence of the 
local spin susceptibility. The degeneracy of the dyz and 
dzx orbitals may play a role here. We point out that 
the essential information should be obtained if the cal- 
culation of the spin density distribution under a uniform 
magnetic field can be made, which we leave for our future 
study. 

It has also been reported^ that, under high pressures, 
a variety of charge-ordered phases appear in this system; 
in particular, a different charge-ordered pattern with the 
antiferromagnetic long-range order has been suggested 
to occur above '--^l GPa at low temperatures. Compet- 
ing phases with similar energies have therefore been sug- 
gested to exist in this system, which should further be 
clarified both experimentally and theoretically. 



IV. DISCUSSION 

Our calculated results presented here contains series of 
theoretical predictions, which should be checked by the 
experimental studies. 

Let us first compare our calculated results with the 
results of the photoemission spectroscopy experiment. * 
It has been reported that, above the transition tempera- 
ture, the coherent quasiparticle peak appears at the bind- 
ing energy of £f, ~ — 0.5 eV (where the Fermi energy 
is set to 0), which we find corresponds well to the t2g 
band below the Fermi energy shown in Fig. 3. The large 
incoherent spectral weight also appears at the binding 
energy of £(, ~ 0.5 — 2.5 eV where a large gap with van- 
ishing density of states appears in our calculated results 
(see Fig. 3). This result suggests the importance of elec- 
tron correlations in this material; as has been pointed 
out in Ref.'*, the situation is similar to the case of V2O3, 
for which the single-particle excitation spectra have been 
studied in detailj^^ The large and broad spectral weight 
appears at Sb — 2.5 — 9 eV, which comes mainly from the 
O 2p orbitals. The position and shape of this spectral 
weight agree well with our calculated density of states 
shown in Fig. 3. 

Below the transition temperature, the coherent peak 
near the Fermi energy has been reported to shift largely 



V. SUMMARY 

We have carried out the ab initio electronic structure 
calculation of hoUandite vanadate K2V8O16 to clarify its 
basic electronic states in the weak correlation limit. We 
have shown the following: 

(i) The states near the Fermi level consist predominantly 
of the three t2g orbitals and the hybridization with the 
oxygen 2p orbitals is small. The strong-coupling model 
based on the three t2g orbitals may therefore be justified. 

(ii) The dyz and dzx orbitals are exactly degenerate and 
are lifted from the dxy orbital. The doubly-degenerate 
atomic energy levels may play a role in the observed elec- 
tronic properties of the system. 

(iii) The calculated band structure and Fermi surface in- 
dicate that the system is not purely one-dimensional but 
the coupling between the VO double chains is important. 

(iv) No nesting features in the Fermi surfaces that con- 
tribute to the instability of the doubling of the unit cell 
along the c axis are seen, suggesting the observed lattice 
instability to be the strong-coupling origin. 

(v) The presence of strong electron correlations is sug- 
gested in the observed valence-band photoemission spec- 
tra compared with the calculated density of states. 

We hope that the present results obtained in the weak 
correlation limit will be combined with the theory based 
on the strongly-correlated electron models in near future. 
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to help us understand the nature of the charge, orbital, 
and spin degrees of freedom of this intriguing material. 
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